£ 



o> 



X 

in: 



(N 
> 

O 
On 
(N 



X 



DESY-08-093 
ZEUS-pub-08-006 
July 2008 



0\ 
O 
O 

^ Production of excited charm and 



charm-strange mesons at HERA 



ZEUS Collaboration 



Abstract 



The production of excited charm, .D 1 (2420)° and D2(2A60)°, and charm-strange, 
D s i(2536) ± , mesons in ep collisions was measured with the ZEUS detector at 
HERA using an integrated luminosity of 126 pb" 1 . Masses, widths and helicity 



parameters were determined. The measured yields were converted to the rates 
of c quarks hadronising as a given excited charm meson and to the ratios of 



the dominant Z^2(2460) and D s i(2536) ± branching fractions. A search for the 
radially excited charm meson, D*'(2640) ± , was also performed. The results are 
compared with those measured previously and with theoretical expectations. 
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1 Introduction 

Heavy-quark spectroscopy has recently undergone a renaissance with the discovery of 
several new states [1]. The properties of these states challenge the theoretical description 
of heavy-quark resonances. Therefore, further measurements of excited charm and charm- 
strange mesons are important. 

The lowest-mass states of the cq (cq) system (q = u, d, s) with spin zero (P mesons) and 
spin one (P* mesons) and with orbital angular momentum L = are well established [1]. 
A singlet and a triplet of states with L = 1 are expected. These P-wave (L — 1) 
mesons can decay to charm mesons with L = by emitting a pion or a kaon. Heavy 
Quark Effective Theory [2] (HQET) predicts that, in the heavy-quark limit (ttiq— >oo), the 
properties of the P-wave mesons are determined mainly by the total angular momentum 
of the light quark, j = L + s, where s denotes the spin of the light quark. Consequently, 
the four states are grouped in two doublets with j = 3/2 or 1/2. Only P-wave decays 
are allowed for the members of the j = 3/2 doublet; therefore they are supposed to be 
narrow. On the other hand, the members of the j = 1/2 doublet decay through 5-wave 
only and therefore are expected to be broader [3]. Due to the finite charm quark mass a 
separation of the two doublets is only an approximation and amplitudes of two observable 
states with J p = 1 + can be mixtures of P- and S'-wave amplitudes. Here J and P are 
the total angular momentum and parity of the cq system. 

Two pairs (neutral and charged) of narrow non-strange excited charm mesons, P!(2420) 0,± 
and P2(2460) 0,± , and a pair of narrow charm-strange excited mesons, P s i(2536) ± and 
P S 2(2573) ± , were observed and tentatively identified as the members of the j = 3/2 
doublets with J p = 1 + and 2 + , respectively [1]. Recently, the HQET expectations were 
supported by the first measurements of the broad non-strange excited charm mesons: 
neutral and charged P*(2400)°' ± with J p = 0+ [4,5], and Pi(2430)° with J p = 1+ [4]. 
The predicted broad non-strange charged excited charm meson with J p = 1 + has not 
yet been observed. The recent discovery of two additional charm-strange excited mesons, 
P s * (2317) ± with J p = 0+ and P sl (2460) ± with J p = 1+ reported initially by BABAR [6] 
and CLEO [7], respectively, revealed their surprisingly small masses and narrow widths [1]. 
The small mass values forbid their decay into D^K final states. 

In addition to the orbital excitations, radially excited charm mesons D'(J P = 0~) and 
D*'(J P = 1~) were predicted with masses of about 2.6 GeV and dominant decay modes 
to Dim and P*7T7t, respectively [8,9]. An observation of a narrow resonance in the final 
state P* ± 7r + 7r~ at 2637 MeV was reported and interpreted as the radially excited P* /=t 
meson by DELPHI [10]. However, OPAL found no evidence for this narrow resonance in 
an analogous search [11]. 

Production of non-excited charm and charm-strange hadrons was extensively studied at 



HERA [12,13]. The large charm production cross section at HERA also provides a means 
to study excited charm and charm-strange mesons produced in ep collisions. The first 
such study is reported in this paper. It is restricted to decays, for which significant signals 
were identified: 



Z>i(2420)° - 


-+ D* + n-, 


L>;(2460)° - 


■+ D* + 7r~ , D + n- . 


L> s i(2536) + - 


■+ D* + K° S ,D*°K 



The corresponding antiparticle decays were also measured 1 . A search for the radially 
excited charm meson, D*'(26A0) + , in the D* + 7i + tt~ final state was also performed. 

2 Experimental set-up 

The analysis was performed using data taken with the ZEUS detector from 1995 to 2000. 
In this period, HERA collided electrons or positrons 2 with energy E e = 27.5 GeV and 
protons with energy E p = 820 GeV (1995-1997) or E p = 920 GeV (1998-2000). The data 
used in this analysis correspond to an integrated luminosity of 126.5 ± 2.4 pb~ . 

A detailed description of the ZEUS detector can be found elsewhere [14]. A brief outline 
of the components most relevant to this analysis is given below. 

Charged particles were tracked in the central tracking detector (CTD) [15], which operated 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consisted of 72 cylindrical drift chamber layers, organized in nine superlayers covering the 
polar-angle 3 region 15° < 9 < 164°. The transverse-momentum resolution for full-length 
tracks was a(p T )/p T = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. To estimate the 
energy loss per unit length, dE/dx, of charged particles in the CTD [16,17], the truncated 
mean of the anode-wire pulse heights was calculated, which removes the lowest 10% and 
at least the highest 30% depending on the number of saturated hits. The measured dE/dx 
values were corrected for a number of effects [18] and normalised such that the corrected 
value was one for a minimum ionising particle. The resolution of the dE/dx measurement 
for full-length tracks was about 9%. 

The high-resolution uranium-scintillator calorimeter (CAL) [19] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part was 



1 Hereafter, charge conjugation is implied. 

2 From now on, the word "electron" is used as a generic term for electrons and positrons. 

3 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 



subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). 
The smallest subdivision of the calorimeter was called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, were a(E)/E = 0.18/ yE for electrons and 
a(E)/E = 0.35/VE for hadrons, with E in GeV. 

The luminosity was determined from the rate of the bremsstrahlung process ep — > e'yp, 
where the photon was measured with a lead-scintillator calorimeter [20] located at Z = 
-107 m. 



3 Event simulation 

Monte Carlo (MC) samples of charm and beauty events were produced with the Pythia 
6.156 [21] and Rapgap 2.0818 [22] event generators. The Rapgap MC used Hera- 
cles 4.6.1 [23] in order to incorporate first-order electroweak corrections. The generation 
included direct photon processes, in which the photon couples directly to a parton in the 
proton, and resolved photon processes, where the photon acts as a source of partons, one 
of which participates in the hard scattering process. The CTEQ5L [24] and GRV LO [25] 
parametrisations were used for the proton and photon structure functions, respectively. 
The charm and bottom quark masses were set to 1.5 GeV and 4.75 GeV, respectively. 
Events for all processes were generated in proportion to the MC cross sections. The Lund 
string model [26] as implemented in Jetset [21] was used for hadronisation in Pythia 
and Rapgap. The Bowler modification [27] of the Lund symmetric fragmentation func- 
tion [28] was used for the charm and bottom quark fragmentation. To generate D*' + 
mesons, which are not present in the Jetset particle table, the mass of a charged charm 
meson in the table was set to 2.637 GeV, its width was set to 15MeV and the decay 
channel was set to D* + 7t + tt~ [10]. 

The Pythia and Rapgap generators were tuned to describe the photoproduction and the 
deep inelastic scattering (DIS) regimes, respectively. Consequently, the Pythia events, 
generated with Q 2 < 0.6 GeV 2 , were combined with the Rapgap events, generated with 
Q 2 > 0.6 GeV 2 , where Q 2 is the exchanged-photon virtuality. Diffractive events, charac- 
terised by a large rapidity gap between the proton at high rapidities and the centrally- 
produced hadronic system, were generated using the Rapgap generator in the diffractive 
mode and combined with the non-diffractive MC sample. The contribution of diffractive 
events was estimated by fitting the rj max distribution 4 of the data with a linear combination 
of the non-diffractive and diffractive MC samples. 



4 The quantity ?7max is denned as the pseudorapidity of the CAL energy deposit with the lowest polar 
angle and an energy above 400 MeV. 



To ensure a good description of the data, the transverse momenta, pt{D* + , D + , D°), and 
pseudorapidity, r](D* + ,D + ,D°), distributions were reweighted to the data for the com- 
bined Pythia+Rapgap MC sample. The reweighting factors, tuned using a large D* + 
sample (Section 4), were used for D + and D° mesons relying on the MC description of the 
differences between the D* + and D + or D° distributions. The effect of the reweighting on 
the measured values was small; the reweighting uncertainty was included when evaluating 
systematic uncertainties (Section 8). 

The generated events were passed through a full simulation of the detector using Geant 
3.13 [29] and processed with the same reconstruction program as used for the data. 



4 Event selection and reconstruction of lowest-mass 
charm mesons 

Events from both photoproduction [30] and DIS [13] were selected online with a three- 
level trigger [14,31]. The first- and second-level trigger used CAL and CTD data to 
select ep collisions and to reject beam-gas events. At the third level, where the full 
event information was available, the nominal charm-meson trigger branches required the 
presence of a reconstructed D* + , D + or D° candidate. The efficiency of the online charm- 
meson reconstruction, determined relative to the efficiency of the offline reconstruction, 
was above 95%. Events missed by the nominal charm- meson triggers but selected with 
any other trigger branch, dominantly from an inclusive DIS trigger and a photoproduction 
dijet trigger, were also used in this analysis. 

In the offline analysis, only events with \Z vtx \ < 50cm, where Z vtx is the primary vertex 
position determined from the CTD tracks, were used. The D* + , D + and D° mesons 
were reconstructed using tracks measured in the CTD and assigned to the reconstructed 
primary event vertex. To ensure both good track acceptance and good momentum reso- 
lution, each track was required to have a transverse momentum greater than 0.1 GeV and 
to reach at least the third superlayer of the CTD. 

To suppress the combinatorial background, a cut on the ratio pt(D* + ,D + ,D°)/E^ >10 , 
motivated by the hard character of charm fragmentation, was applied. The transverse 
energy, E^ lQ ° , was calculated as Ei i ^>io°(£'iSin^), where the sum runs over all energy 
deposits in the CAL with the polar angle 9 outside a cone of 9 = 10° around the forward 
direction. Moreover, the measured dE/dx values of those tracks that were candidates to 
come from D* + , D + and D° were used. The parametrisations of the dE/dx expectation 
values and the \ 2 probabilities Ik and /„- of the kaon and pion hypotheses, respectively, 
were obtained in the same way as described in previous publications [30,32]. To maximise 



the ratios of the numbers of correctly assigned kaons and pions to the square roots of the 
numbers of background particles, the cuts Ik > 0.03 and l n > 0.01 were applied. 

The measurements were done in the full kinematic range of Q 2 . Events produced in the 
photoproduction regime with Q 2 < 1 GeV 2 contributed 70 — 80 % of the selected D* + , D + 
and D° samples. 

4.1 Reconstruction of D* + mesons 

The D* + mesons were identified using the two decay channels 

D *+ _> D o n + _» (jr*r>+ (1) 

D*+ -► P°tt+ -► (A'"7r + 7r + 7T-)7r+. (2) 

The pion from the D* + decays is referred to as the "soft" pion, tt s , because it is constrained 
to have limited momentum by the small mass difference between the D* + and D° [1]. 

Selected tracks were combined to form D° candidates assuming the decay channels (1) or 
(2). For both cases, D° candidates were formed by calculating the invariant mass M(Kir) 
or M{Kirinr) for combinations having a total charge of zero. The soft pion was required 
to have a charge opposite to that of the particle taken as a kaon and was used to form 
a D* + candidate having mass M(Kinr s ) or M(Kinnnr a ). To reduce the combinatorial 
background, requirements (see Table 1) similar to those used in a previous publication [32] 
were applied. 

The mass difference AM = M(Kinr s ) - M(Ktt) for channel (1) or AM = M(Kinnnr s ) - 
M(K7nnr) for channel (2) was evaluated for all remaining D* + candidates. Figures la 
and lb show the mass-difference distributions for channels (1) and (2), respectively. Peaks 
at the nominal value of M(D* + ) — M(D°) are evident. 

To determine the background under the peaks, wrong-charge combinations were used. For 
both channels (1) and (2), these are defined as combinations with total charge ±2 for the 
D° candidate and total charge ±1 for the D* + candidate. The histograms in Fig. 1 show 
the AM distributions for the wrong-charge combinations, normalised to the distributions 
of D* + candidates with the appropriate charges in the range 0.15 < AM < 0.1685 GeV 
for channel (1) and 0.15 < AM < 0.16 GeV for channel (2). The upper ends of the 
normalisation ranges correspond to the trigger selections of D* + candidates in the two 
decay channels. The multiple counting of a D* + candidate produced by D° candidates 
formed by the same tracks was excluded [32]. 

To improve the signal-to-background ratio, only D* + candidates with 0.144 < AM < 
0.147 GeV for channel (1) and 0.1445 < AM < 0.1465 GeV for channel (2) were kept 



for the excited charm and charm-strange meson studies. After background subtraction, 
signals of 39500±310 D* + mesons in channel (1) and 17300±210 D* + mesons in channel (2) 
were found in the above AM ranges. 

The AM distributions were also fitted to a sum of a modified Gaussian function describing 
the signal and a background function. The modified Gaussian function was defined as 

Gauss mod oc exp[-0.5 • x 1 ^^ - 5 ^], (3) 

where x = |(AM — Mq)/(t\. This functional form described both data and MC signals 
well. The signal position, M , and width, a, as well as the numbers of D* + mesons in 
the signal window were free parameters of the fit. The background function had a form 
A ■ (AM — m n +) B ■ exp[C • (AM — m n +)], where m n + is the pion mass [1] and A, B and 
C were free parameters. The fit yielded mass difference values of 145.46 ± 0.01 MeV for 
channel (1) and 145.45 ± 0.01 MeV for channel (2), in agreement with the PDG value [1]. 
The widths of the signals were 0.59±0.01 MeV and 0.51±0.01 MeV, respectively, reflecting 
the detector resolution. 

4.2 Reconstruction of D + mesons 

The D + mesons were reconstructed from the decay D + — + K~ii + ii + . In each event, two 
tracks with the same charges and px > 0.5 GeV and a third track with opposite charge and 
Pt > 0.7 GeV were combined to form D + candidates. The pion masses were assigned to 
the two tracks with the same charges and the kaon mass was assigned to the third track, 
after which the candidate invariant mass, M(Ktt7t), was calculated. To suppress the 
combinatorial background, a cut of cos9*(K) > —0.75 was imposed, where 8*(K) is the 
angle between the kaon in the Knir rest frame and the Kttti line of flight in the laboratory 
frame. To further suppress the combinatorial background, a cut pt(D + )/E t >10 > 0.25 
was applied. To suppress background from D* + decays, combinations with M(Kim) — 
M(Ktt) < 0.15 GeV were removed. The background from Df — > <p7i + with <p — » K + K~ 
was suppressed by requiring that the invariant mass of any two D + candidate tracks with 
opposite charges was not within ±8 MeV of the nominal (f) mass when the kaon mass was 
assigned to both tracks. Only D + candidates in the kinematic range px{D + ) > 2.8 GeV 
and —1.6 < i](D + ) < 1.6 were kept for further analysis. 

Figure 2a shows the M(K7nr) distribution for the D + candidates after all cuts. Reflections 
from Df and A+ decays to three charged particles were subtracted using the simulated 
reflection shapes normalised to the Df and A+ production rates previously measured by 
ZEUS [30]. A clear signal is seen at the nominal value of the D + mass. To improve the 
signal-to-background ratio, only D + candidates with 1.850 < M(Ktt7t) < 1.890 GeV were 
kept for the excited charm meson studies. The mass distribution was fitted to a sum of 



a modified Gaussian function describing the signal and a linear function describing the 
non-resonant background. The fit yielded a D + mass value 1867.9 ±0.5 MeV in agreement 
with the PDG value [1]. The width of the signal was 12.9±0.5 MeV, reflecting the detector 
resolution. The number of D + mesons yielded by the fit in the above M(Ktttt) range was 
N(D + ) = 20430 ± 510. 



4.3 Reconstruction of D° mesons 

The D° mesons were reconstructed from the decay D° — > K~ir + . In each event, tracks 
with opposite charges and p? > 0.8 GeV were combined in pairs to form D° candidates. 
To suppress the combinatorial background, a cut of \ cos 0*(K)\ < 0.85 was imposed, 
where 8*(K) is the angle between the kaon in the Ktt rest frame and the Kit line of 
flight in the laboratory frame. To further suppress the combinatorial background, a cut 
p T {D Q )/E e > 1Q ° > 0.25 was applied. 

For selected D° candidates, a search was performed for a track that could be the soft 
pion in a D* + — > D°ir^ decay. The soft pion was required to have pr > 0.1 GeV and a 
charge opposite to that of the particle taken as a kaon. The corresponding D° candidate 
was rejected if the mass difference, AM = M(Knir s ) — M(Kir), was below 0.15 GeV. All 
remaining D° candidates were considered "untagged", i.e. not originating from identified 
D* + decays. Only D candidates in the kinematic range pt{D°) > 2.8 GeV and —1.6 < 
T)(D°) < 1.6 were kept for further analysis. 

Figure 2b shows the M(Kir) distribution for untagged D° candidates after all cuts. A 
reflection, produced by D° mesons with the wrong (opposite) kaon and pion mass assign- 
ment, was subtracted using the rejected sample of the D° mesons originating from D* + 
decays [30] . A clear signal is seen at the nominal value of the D° mass. To improve the 
signal-to-background ratio, only D° candidates with 1.845 < M(Kti) < 1.885 GeV were 
kept for the excited charm-strange meson studies. The mass distribution was fitted to 
a sum of a modified Gaussian function describing the signal and a background function. 
Monte Carlo studies showed that the background shape was compatible with being linear 
in the mass range above the signal. For smaller M(Ktt) values, the background shape 
exhibits an exponential enhancement due to contributions from other D° decay modes 
and other D mesons. Therefore the background shape in the fit was described by the form 
[A + B ■ M(Kn)} for M(Ktt) > 1.86 GeV and [A + B ■ M(Kir)] ■ exp{C ■ [M(Kn) - 1.86]} 
for M(Ktt) < 1.86 GeV, where A, B and C were free parameters. The fit yielded the 
D° mass value 1864.9 ± 0.2 MeV in agreement with the PDG value [1]. The width of the 
signal was 17.4 ±0.2 MeV, reflecting the detector resolution. The number of untagged D° 
mesons yielded by the fit in the above M(Ktt) range was N(D^ nt ) = 22110 ± 440. 



5 Study of the excited charm mesons D® and D^ 

5.1 Reconstruction of D®, D*° —*■ D* + n~ decays 

To reconstruct the D®, D^ — > D* + ir~ decays, an excited charm meson candidate was 
formed by combining each selected D* + candidate (Section 4.1) with an additional track, 
assumed to be a pion (ir a ), with a charge opposite to that of the D* + candidate. The addi- 
tional track was required to satisfy the pion dE/dx hypothesis with l n > 0.01 (Section 4). 
To reduce the combinatorial background, the following requirements were applied: 

PtM > 0.15 GeV, p T (D* + ir a )/E^ >10 ° > 0.25, cos 0* (£>*+) < 0.9 

for the D* + decay channel (1), and 

p T (ir a ) > 0.25 GeV, p T (D*+-K a )/E^ >w ° > 0.30, cos 0* (£>*+) < 0.8 

for the D* + decay channel channel (2). The decay angle 6*(D* + ) is the angle between 
the D* + in the D* + ir a rest frame and the D* + ir a line of flight in the laboratory frame. A 
cut ?7(7r a ) < 1.1 was applied to exclude the region of large track density in the forward 
(proton) direction. 

For each excited charm meson candidate, the "extended" mass difference, AM ext = 
M(Kirir s iT a ) - M(Kim s ) or AM cxt = M(Kiiinm s ir a ) - M(Kirinm s ), was calculated. The 
invariant mass of the D* + ir a system was calculated as M(D* + ir a ) = AM cxt + M(D* + )p DG , 
where M(D* + )pdq is the nominal D* + mass [1]. The resolution in M(D* + ir a ) around the 
nominal masses of the D® and D^ mesons [1] was estimated from MC simulations to be 
5.6 MeV. 

Figure 3a shows the M(D* + ir a ) distribution for D* + meson candidates reconstructed 
in both decay channels (1) and (2). A clear enhancement is seen in the range 2.4 < 
M(D* + ir a ) < 2.5 GeV, where contributions from Di(2420)° and ^2(2460)° mesons are ex- 
pected. The wide _Di(2430)° meson, which is also expected to contribute to the M(D* + ir a ) 
distribution, is not distinguishable from background due to its large width (384" I 1 1 .J?J ± 
74 MeV [1]). No enhancement is seen in the M(D* + n a ) distribution for wrong charge 
combinations (histogram) formed by combining a D* + candidate and ir a with the same 
charges. The wrong charge distribution lies generally below the distribution for the combi- 
nations with the appropriate charges, in agreement with MC predictions; this is expected 
near threshold since, due to charge conservation, the invariant mass distribution for ran- 
dom track combinations with total charge ±2 should lie below that for track combinations 
with total charge zero. 



5.2 Reconstruction of D*° — > D + 7r decays 

To reconstruct the D^ — > D + tt~ decays, an excited charm meson candidate was formed 
by combining each selected D + candidate (Section 4.2) with an additional track, assumed 
to be a pion (vr a ), with a charge opposite to that of the D + candidate. The additional 
track was required to satisfy the pion dE/dx hypothesis with l n > 0.01 (Section 4). To 
reduce the combinatorial background, the following requirements were applied: 

ri(n a ) < 1.1, p T (n a ) > 0.30 GeV, p T (D+7r a )/ ' E e T >w > 0.35, cos6*(D + ) < 0.8, 

where 6*(D + ) is the angle between the D + in the D + ir a rest frame and the D + n a line of 
flight in the laboratory frame. 

For each excited charm meson candidate, the extended mass difference, AM ext = 
M(i^7r7T7r a ) — M(KTTir), was calculated. The invariant mass of the D + n a system was 
calculated as M(D + ir a ) = AM cxt + M(D + ) PDG , where M(D + ) PDG is the nominal D + 
mass [1]. The resolution in M(D + 7i a ) around the nominal mass of the D^ meson [1] was 
estimated from MC simulations to be 7.3 MeV. 

Figure 3b shows the M(D + 7r a ) distribution for the selected excited charm meson can- 
didates. A small excess is seen around the nominal mass of the D^ meson. The wide 
-Dq(2400)° meson, which is also expected to contribute to the M(D + ir a ) distribution, is 
not distinguishable from background due to its large width (261 ± 50 MeV [1]). As ex- 
pected from parity and angular momentum conservation for a 1 + state, no indication 
of the D\ decay to D + ix~ is seen. Feed-downs from the D® and D^ mesons decaying 
to D* + ix~ with a consequent D* + decay to a D + and undetected neutrals, predicted by 
MC at M(D + 7T a ) ~ 2.3 GeV, are not seen, probably due to the large combinatorial back- 
ground. No signal is seen in the M(D + ir a ) distribution for wrong charge combinations 
(histogram) formed by combining a D + candidate and a 7r a with the same charges. 



5.3 Mass, width and helicity parameters 

To distinguish the D\ (1 + state from j = 3/2 doublet) and Df (2 + state from j = 3/2 
doublet) mesons from each other and from the wide _Di(2430)° (1 + state from j = 1/2 
doublet) meson, the helicity angular distribution was used. The helicity angle (a) is 
defined as the angle between the n a and ir s momenta in the D* + rest frame. The helicity 
angular distribution can be parametrised as 

dN , 2 

oc 1 + h cos a, (4) 

a cos a 

where h is the helicity parameter. HQET predicts h — 3 (h — 0) for the 1 + state from 
the j = 3/2 (j = 1/2) doublet, and h = —1 for the 2 + state from the j = 3/2 doublet. 



Figure 4 shows the M(D* + ir a ) distribution in four helicity intervals. The Dj-meson contri- 
bution is increasing with | cos(a)| and dominates the excess in the M(D* + 7r a ) distribution 
for | cos(a)| > 0.75. The dependence of the D^-meson contribution on the helicity angle 
is less pronounced; it is consistent with the expected slow decrease with | cos(a)|. 

To extract the D® and Dl^ yields and properties, a minimal \ 2 fit was performed using si- 
multaneously the M(D + ir a ) distribution (Fig. 3b) and the M(D* + ir a ) distributions in four 
helicity intervals (Fig. 4). Each of the D? -> D* + n-, Df -» D* + tt~ and Df -► D+n- 
signals was represented in the fit by a relativistic D-wave Breit-Wigner function (see 
Appendix) convoluted with a Gaussian resolution function with a width fixed to the cor- 
responding MC prediction. The dependence of the detector acceptance and resolution on 
the M(D* + ir a ) or M(D + ir a ) was obtained from MC and corrected for in the fit function. 
Equation (4) was used to describe the helicity distributions. The acceptance dependence 
on the helicity angle, found from MC to be very weak, was corrected for in the fit func- 
tion. Yields of all three signals, the D® and Df masses, and the D\ width and helicity 
parameters were free parameters of the fit. Since the data were not able to constrain 
reliably the Df width and helicity parameter, the Df width was fixed to the recently 
updated world average value of 43 ± 4MeV [1] and the HQET prediction, h(Df) = — 1, 
was used for the helicity parameter. 

To describe backgrounds in the M(D* + 7r a ) and M(D + 7T a ) distributions, a functional form 
with three shape parameters x A exp(— Bx + Cx 2 ), where x = AM ext — m n +, was used. 
It was checked that such a functional form describes the wrong charge distributions well. 
The yields and shape parameters of the M(D* + n a ) and M(D + n a ) background functions 
were independent free parameters of the fit. Since neither data nor MC demonstrated a 
sizeable background dependence on the helicity angle, the same background function was 
used for the M(D* + -ir a ) distributions in the four helicity intervals. 

The expected feed-downs from D®,Df — > D* + n~ —* D + n~+ neutrals (Section 5.2) were 
included in the M (D + 7r a ) fit function; the effect on the fit results was small. Contributions 
from the wide L>i(2430) and £>*(2400)° states were added to the M(D* + ir a ) and M(D + 7r a ) 
fit, respectively. Their shapes were described with a relativistic S-wave Breit-Wigner 
function (see Appendix) convoluted with a Gaussian resolution function with widths fixed 
to the MC prediction. The masses and widths of the wide excited charm mesons were 
set to the world-average values [1]. The _Di(2430)° yield was set to that of the narrow 
-D!(2420)° meson since both have the same quantum numbers. The -Dq(2400)° yield was 
set to 1.7 times the Df — > D + tt~ yield as observed by the FOCUS collaboration [5]. 
The yield measured by FOCUS covers both a direct signal from the -Dq(2400)° and a 
feed-down from the _Di(2430)°, decaying to D* + tt~ with a consequent D* + decay to a D + 
and undetected neutrals [5]. 

The results of the simultaneous fit including all contributions are shown in Figs. 3-4. The 
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fit with 15 free parameters described well the M(D + ir a ) distribution and the M(D* + 7i a ) 
distributions in four helicity intervals with a \ 2 °f 913 for 925 degrees of freedom. The 
numbers of reconstructed excited charm mesons and values of all free background param- 
eters yielded by the fit are summarised in Table 2. 

The differences between the D\ and D^ masses and M(D* + )pdg were 

M{Dl) - M(D*+) PDG = 410.2 ± 2.1(stat.) ± 0.9(syst.) MeV, 

M(Df) - M(D* + ) PBG = 458.8 ± 3.7(stat.)+JJ(syst.) MeV, 
and, hence, the masses of the D® and D^ were 

M(Dl) = 2420.5 ± 2.1(stat.) ± 0.9(syst.) ± 0.2(PDG) MeV, 

M(Df) = 2469.1 ± 3.7(stat.)ii|(syst.) ± 0.2(PDG) MeV. 

The first uncertainties are statistical, the second are systematic (Section 8) and the third 
are due to the uncertainty of the M(D* + ) PDG value. Small errors due to the uncertainty 
of the M(D* + )pdq — M(D + )pdq value were included in the systematic uncertainties. The 
measured D\ and D^ masses are in fair agreement with the world average values [1]. The 
D\ width yielded by the fit is 

r(D?) = 53.2 ± 7.2(stat.)+^(syst.) MeV 

which is above the world average value 20.4 ± 1.7 MeV [1]. The observed difference can 
be a consequence of differing production environments. The D® width can have a sizeable 
contribution from the broad S- wave decay even if the 5"- wave admixture is small [33,34]. 
A larger S-wave admixture at ZEUS with respect to that in measurements with restricted 
phase space, which can suppress production of the broad state, could explain why the 
measured D® width is larger than the world average value. 

The D® helicity parameter was 

h{D\) = 5.9l?;?(stat.)i?:o(syst.). 

This is inconsistent with the prediction for a pure S'-wave decay of the 1 + state, h = 0. 
It is consistent with the prediction for a pure .D-wave decay, h = 3. 

In the general case of D- and S'-wave mixing, the helicity angular distribution form of the 
1 + state is: 

ocr + (l -r)(l + 3cos 2 a)/2+ ^/2r(l - r) cos 0(1 -3 cos 2 a), (5) 



dcosa 
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where r = Ts/{Ts H-To), Ts/d is the S-/D-wave partial width and <p is the relative phase 
between the two amplitudes. Using Eqs. (4) and (5), cos</> can be expressed in terms of 
r and the measured value of the helicity parameter, h: 

(3-h)/(3 + h)-r 

COS0= = . (6) 

2y/2r(l-r) 

Figure 5 compares with previous measurements the range restricted by the measured 
h(D\) value and its uncertainties in a plot of cos</> versus r. The ZEUS range has a 
marginal overlap with that restricted by the CLEO measurement of h(Df) = 2.74^Q'g3 [35]. 
BELLE performed a three-angle analysis and measured both the cos0 and r values [4]. 
The BELLE measurement, which suggested a very small admixture of S'-wave to the 
Di(2420) — > D* + n~ decay and almost zero phase between two amplitudes, is outside the 
ZEUS range; the difference between the two measurements, evaluated with Eq. (6), is ~ 2 
standard deviations. 



5.4 Fragmentation and branching fractions 

The numbers of reconstructed D®, D^ — * D* + ir~ and D^ — > D + n~ decays were divided 
by the numbers of reconstructed D* + and D + mesons, yielding the rates of D* + and 
D + mesons originating from the D® and D^ decays. To correct the measured rates 
for detector effects, the relative acceptances were calculated using the MC simulation as 
ratios of acceptances for the D®, D^ —* D* + tt~ and D^ —* D + -k~ states to the inclusive 
D* + and D + acceptances, respectively. The acceptance of the requirement l v > 0.01 
for the additional track was calculated with data using identified pions from D* + decays 
(Section 4.1), to be (98.9 ± 0.1)%; only pions in the kinematic range of the additional 
pion selection were used. 

Charm production at HERA is larger than beauty production by two orders of magni- 
tude. The small 6-quark relative contributions, predicted by the MC simulation using 
branching fractions of 6-quark decays to the charm hadrons measured at LEP [36-39] 5 , 
were subtracted when calculating the relative acceptances; the subtraction changed the 
relative acceptances by less than 1.5% of their values. The relative acceptances were 52% 
for the D®, D^ —>■ D* + ix~ and 47% for D^ — > D + ir~ in the kinematic ranges described 
in Section 4. 

The fractions, JF, of D* + mesons originating from D® and D^ decays were calculated in 
the kinematic range \r](D* + )\ < 1.6 and pt{D* + ) > 1.35 GeV for the D* + decay channel 



5 The published branching fractions of the 6-quark decays were recalculated using updated values [1] of 
the relevant charm-hadron decay branching fractions. 
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(1), combined with channel (2) for p T (D* + ) > 2.8 GeV: 

7 : d°->d* + *-/d*+ = 10-4 ± 1.2(stat.)l?|(syst.) %, 
^Df^D*+n-/D*+ = 3.0 ± 0.6(stat.) ± 0.2(syst.) %. 

The fraction of D + mesons originating from Df decays, calculated in the kinematic range 
p T (D+) > 2.8 GeV and \rj(D+)\ < 1.6 is 

Fd^d+*-/d+ = 7.3 ± 1.7(stat.)lS:i(syst.) %. 

The fractions measured in the restricted pt{D* + , D + ) and rj(D* + , D + ) kinematic ranges 
were extrapolated to the fractions in the full kinematic phase space using the Bowler 
modification [27] of the Lund symmetric fragmentation function [28] as implemented in 
Pythia [21]. Applying the estimated extrapolation factors, ~ 1.1 for J 7 d°,d*°^d*+k- m*+ 
and ~ 1.2 for F d ?-*d+-k-/d+> gives 

Ftf>->D* + *-/D* + = H-6 ± 1.3(stat.)+^(syst.) %, 

FDf^D*+n-/D«+ = 3-3 ± 0.6(stat.) ± 0.2(syst.) %, 

F%?->d + *-/d + = 8-6 ± 2.0(stat.)ili(syst.) %. 



In the full kinematic phase space, the extrapolated fractions of D* + originating from D® 
and Df and of D + originating from Df can be expressed as 

-r-extr _ J\ C ~* D\) K? 

•' nl) ,. n* + „- /n*+ — „ , „.,. i n ' O l 



' D1^D*+tt-/D*+ ~ f / n *+\ - D D°^D*+n 



f(c~ 


■+ D*+) 


f(c~ 


- Df) 


f(c- 


■* D*+) 


f(c~ 


♦ Df) 



<£-extr J y- ~_z/ 10 

- r D*°^D*+7r-/D*+ _: 777^7^+^ ' ,D ^ ^£)* + 'r-' 

^D|°-»D+7r-/D+ - y7T77£J+y ' D ^°-»-D+7r-' 

where the fragmentation fractions /(c — > D®), f(c — > Df), f(c — > £)* + ) and /(c — » D + ) 
are the rates of c quarks hadronising as a given charm meson, and B D o_ >D , +7T - , £> D »o^ D , +7r - 
and £> D *o^ D+7r - are the corresponding branching fractions. 

These expressions provide a means to calculate the fragmentation fractions f(c — > _D°) 
and /(c — > -D^ ), and the ratio of the two branching fractions for the Df meson: 

-rextr 
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— D*+7T- 
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The /(c — > -D* + ) and /(c — > £) + ) values, previously measured by ZEUS [30], were recal- 
culated with the updated PDG values of the branching fractions [1] to be 

f(c -► D*+) = 20.4±0.9(stat.)lS:r(syst.)l?:I(br.) %, 

/(c _+ £>+) = 21.7 ± 1.4(stat.)l^(syst.)i 1 1 i(br.) %, 
where the third uncertainties are due to the branching-fraction uncertainties. This yields 

Df^D+n- = 2 g ± .8(stat.)io6(syst.) 
I3d* ^d*+7t- 

in agreement with the world average value of 2.3 ± 0.6 [1]. Theoretical models [34,40,41] 
predict the ratio to be in the range from 1.5 to 3. 

Assuming isospin conservation, for which 

&D°->D*+ir- = 2/3, &D*°-*D*+ir- + #D*°-*£>+*- = 2 /3, 

yields f(c — > D®) and f(c — » Df) (Table 3). In order to check fragmentation universal- 
ity for the excited charm mesons, the measured fragmentation fractions are compared 
and found to be consistent with those obtained in e + e~ annihilations. The measured 
f(c — ► D®) and f(c — » Df) values are above the predictions of the thermodynamical 
model [42] (Table 3). The sum of the two fragmentation fractions, 

f(c -> D?) + f(c - Df) = 7.3 ± 0.8(stat.)+^(syst.) %, 

agrees with the prediction of the tunnelling model of 8.5% [43]. The predictions of both 
models are based on fits to the production rates of light-flavoured hadrons at LEP. 

The ratio 

f(c -> D?)//(c -> Df) = 0.93 ± 0.20(stat.) ± 0.16(syst.) 

is consistent with the simple spin-counting prediction of 3/5. Both thermodynamical and 
tunnelling models suggest the ratio should exceed the spin-counting prediction due to the 
difference between the D\ and Df masses. 
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6 Study of the excited charm- strange meson D^ 

6.1 Reconstruction of D^ — ► D* + Kg decays 

The Kg mesons were reconstructed in their charged-decay mode, Kg — * 7r + 7r~, for those 
events containing a D* + candidate. To identify Kg candidates, displaced secondary ver- 
tices reconstructed from pairs of oppositely charged tracks [44] were used. The identifica- 
tion efficiency degraded for the displaced secondary vertices close to the primary vertex. 
Therefore, additional secondary vertices were formed from pairs of oppositely charged 
tracks that were not assigned to one of the displaced secondary vertices. This was done 
by calculating the intersection points of the two tracks in the XY plane and requiring 
|AZ| < 3 cm between the two tracks at the intersection point. To reduce the combinato- 
rial background originating from tracks from the primary vertex, the additional secondary 
vertices with distances between the primary and secondary vertices in the XY plane of 
less than 0.5 cm were removed. 

To reduce the combinatorial background, it was required that pt > 0.15 GeV for each 
track from any Kg candidate, cosa XY > 0.97 and cosa^ z > 0.85, where a XY and a^ z are 
the projected angles in the XY and 0Z planes, respectively, between the i^-candidate 
momentum and the line joining the primary to the secondary vertex. Figure 6 shows 
the invariant-mass, M(ir + 7T~), distribution for all remaining Kg candidates. Only Kg 
candidates with 0.480 < M{n + n~) < 0.515 GeV were kept for the reconstruction of 
excited charm-strange mesons. The mass distribution was fitted to a sum of a modified 
Gaussian function describing the signal and a linear function describing the non-resonant 
background. The fit yielded the Kg mass value 497.8 ± 0.1 MeV, in agreement with the 
PDG value [1]. The width of the signal was 4.1±0.1 MeV reflecting the detector resolution. 
The number of reconstructed Kg mesons in the range 0.480 < M(ir + 7T~) < 0.515 GeV 
yielded by the fit was N(K° S ) = 8540 ± 120. 

To reconstruct the D^ x — > D* + Kg decays, a .D^-meson candidate was formed by com- 
bining each selected D* + candidate (Section 4.1) with the Kg candidates reconstructed 
in the same event. For each Df x candidate, the extended mass difference, AM ext = 
M(K7Tir s 7i + 7i-) - M(KinT s ) - M(7r + 7r") or AAf cxt = M(KimTT7T s 7T + 7T-) -M(KTTirinT s ) - 
M(7r + 7T~), was calculated. The invariant mass of the D* + Kg system was calculated as 
M(D* + K° S ) = AM cxt + M(D*+) PDG + M(K°) PDG , where M(K° S ) PDG is the nominal K° s 
mass [1]. The resolution in M(D* + Kg) around the nominal mass of the Df t [1] was 
estimated from MC simulations to be 2.2 MeV. 

Figure 7a shows the M(D* + Kg) distribution for D* + meson candidates reconstructed in 
both decay channels. A clear signal is seen at the nominal value of M(D^ X ). 



15 



6.2 Reconstruction of D^ — >■ D*°K + decays 

Monte Carlo studies show that a signal from the DJ^ —> D*°K + decay, with a consequent 
D*° decay to a D° and undetected neutrals, should be seen in the M(D°K + ) distribution 
with an average negative shift of 142.4±0.2 MeV with respect to the nominal D^ mass [1], 
and that the shape of the signal can be reasonably well described by the modified Gaussian 
function (Eq. 3) with a width of 3.1 MeV. 

To reconstruct the D^ x — > D*°K + decays, an excited charm-strange meson candidate was 
formed by combining each selected untagged D° candidate (Section 4.3) with an additional 
track, assumed to be a kaon (K a ), with a charge opposite to that of the particle taken as 
a kaon to form the D° candidate. The additional track was required to satisfy the kaon 
dE/dx hypothesis with l K > 0.03 (Section 4). To reduce the combinatorial background, 
the following requirements were applied: 

r]{K a ) < 1.1, p T (K a ) > 0.60 GeV, p T (D°K a )/E°> 10 > 0.35, cos6*(D°) < 0.8, 

where 6*(D°) is the angle between the D° in the D°K a rest frame and the D°K a line of 
flight in the laboratory frame. 

For each excited charm-strange meson candidate, the extended mass difference, AM ext = 
M(KnK a ) — M(Ktt) was calculated. The invariant mass of the D°K a system was calcu- 
lated as M(D°K a ) = AM ext + M(D°) PDG , where M(D°) PDG is the nominal D° mass [1]. 

Figure 7b shows the M(D°K a ) distribution for the selected excited charm-strange meson 
candidates. A signal is seen at the expected position of the feed-down from the Df x — > 
D*°K + decay. No signal from the known decay D s2 (2573) + — > D°K + [1] was observed, 
probably due to the large combinatorial background. 

6.3 Mass, width and helicity parameters 

The M(D* + Kg) distribution in four helicity intervals is shown in Fig. 8, with the helicity 
angle (a) defined as the angle between the K$ and ir s momenta in the D* + rest frame. 
The D^ x signal decreases with | cos(a)|. 

To extract the D^ yields and properties, an unbinned likelihood fit was performed using 
simultaneously values of M(D°K a ), M(D* + K$), and cos(a) for D* + Kg combinations. 
The observed narrow signals in the M(D* + Kg) and M(D°K a ) distributions were de- 
scribed in the fit by a Gaussian function and a modified Gaussian function, respectively. 
Equation (4) was used to describe the helicity distribution. The acceptance dependence 
on the helicity angle, found from MC to be very weak, was corrected for in the fit func- 
tion. The average shift of the signal in the M(D°K a ) distribution with respect to the 
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mass of D^ meson was fixed to the MC prediction (Section 6.2). Yields and widths of 
both signals, the D^ mass and the D^ helicity parameter were free parameters of the fit. 

To describe the background in the M(D* + K^) distribution, a function x A , where x = 
AM ext , was used. The background description for the M(D°K a ) distribution required 
a functional form with two shape parameters x A exp(— Bx), where x = AM cxt — m^+ 
and rriK+ is the kaon mass [1]. The shape parameters of the M(D* + Kg) and M(D°K a ) 
background functions were independent free parameters of the fit. Since neither data nor 
MC demonstrated a sizeable background dependence on the helicity angle, the background 
function for D* + K s combinations was assumed to be helicity independent. The numbers 
of reconstructed D^ mesons and values of all free background parameters yielded by the 
fit are summarised in Table 4. 

The widths of both signals yielded by the fit agree with the MC predictions for the 
corresponding resolutions. Thus the value of the natural D^ width cannot be extracted. 
The difference between the D^ mass and M(D* + )pdq was 

M(£>+) - M(D*+) PDG = 525.30to;^(stat.) ± 0.10(syst.) MeV, 

and, hence, the mass of the Df ± was 

M(D+) = 2535.57l°;^(stat.) ± 0.10(syst.) ± 0.17(PDG) MeV. 

The first uncertainty is statistical, the second is systematic (Section 8) and the third is 
due to the uncertainty of the M(.D* + )pdg value. Small errors due to the uncertainties 
of the M(D* + ) PDG — M(_D°) PDG and M(K^) PDG values were included in the systematic 
uncertainty. The measured D^ mass is in good agreement with the world average value [1]. 

The D^ helicity parameter was 

h(Dt ± ) = -0.74^(stat.)+^(syst.). 

The measured h value is inconsistent with the prediction for a pure D-wave decay of 
the 1 + state, h = 3, and is barely consistent with the prediction for a pure S'-wave 
decay, h = 0. Figure 9 shows a range, restricted by the measured h^D^) value and 
its uncertainties, in a plot of cos0 versus r = r s /(r s + T D ) (Eq. 6). The measurement 
suggests a significant contribution of both D- and S'-wave amplitudes to the _D s i(2536) + — > 
D* + K® decay. The ZEUS range agrees with that restricted by the CLEO measurement 
of h(Df 1 ) = — 0.23^Q32 [45] and with the BELLE three-angle measurement of both cos0 
and r values [46]. 

6.4 Fragmentation and branching fractions 

The numbers of reconstructed D^ — > D* + Kg and D^ — > D*°K + decays were divided by 
the numbers of reconstructed D* + and untagged D° mesons, respectively, yielding rates of 
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D* + and untagged D° mesons originating from D^ decays. To correct the measured rates 
for detector effects, the relative acceptances were calculated using the MC simulation as 
ratios of acceptances for the D^ — ► D* + K$ and D^ — > D*°K + states to the inclusive 
D* + and untagged-/} acceptances, respectively. The untagged-/} acceptance included 
subtraction of a small contamination to N(D® nt ) from unidentified D* + mesons. The 
acceptance of the requirement Ik > 0.03 for the additional track was calculated with data 
using identified kaons from D* + decays (Section 4.1), to be (95.3 ±0.2)%; only the kaons 
from the kinematic range of the additional kaon selection were used. Subtraction of the 
small 6-quark contribution changed the relative acceptances by less than 2.2% of their 
values. The relative acceptances were 38% for D+ x — > D* + K° s and 48% for D^ — > D*°K + 
in the kinematic ranges described in Section 4. 

The fraction, JF, of D* + mesons originating from D^ x decays, corrected to the fraction 
of K° mesons decaying as Kg (50%) and to the branching fraction of the K s decay 
into 7r + 7r~ (69.20% [1]), was calculated in the kinematic range |?7(.D* + )| < 1.6 and 
Pt{D* + ) > 1.35 GeV for the D* + decay channel (1), combined with channel (2) for 
p T (D* + ) > 2.8 GeV: 



Fd+->d*+k°/d*+ = 1.35±0.18(stat.) ± 0.03(syst.) 

si / 



/o. 



The fraction of untagged D° mesons originating from D^ decays, calculated in the kine- 
matic range pt(D°) > 2.8 GeV and |?7(-D )| < 1.6 is 

Fd+^d^k+zd^ = 1-28 ± 0.26(stat.) ± 0.07(syst.) %. 

The fractions measured in the restricted pt{D* + , D°) and r)(D* + ,D°) kinematic ranges 
were extrapolated to the fractions in the full kinematic phase space (Section 5.4). Ap- 
plying the estimated extrapolation factors, ~ 1.2 for J- ' D o _ > d*+a'°/d*+ an d ~ 1-5 for 
•^D+^o x+/D o ntag , gives 

^d+^+ko/d^ = L67 ± 0-22(stat.) ± 0.07(syst.) %, 

si / 

^d+^d^ k+ /do = L93 ± 0.40(stat.)i°;^( S yst.) %. 

si / untag 

In the full kinematic phase space, the extrapolated fractions of D* + and untagged D° 
mesons originating from Df ± can be expressed as 

•r-cxtr _ j( c ~ i> D s l) n 



D+ l ->D*+K°/D*+ ff c _> £)*H 



j(c — ► D slJ 



'si 



-rextr J \ sij to 



where the fragmentation fractions f(c — > Df x ), f(c — ► D* + ) and /(c — > -Du ntag ) are the 
rates of c quarks hadronising as a given charm meson, and B D + _> D *+ K o and B D + ^ D *o K + 

si si 

are the corresponding branching fractions. 

These expressions provide a means to calculate the fragmentation fraction f(c—> Df x ) 
and the ratio of the two D^[ branching fractions: 



/(c-i>i; 



T extr . f(r — ► D*+"\ -I- 7rextr j/ t-)0 \ 



^n+^n*+!vo + £>? 



>D+-^D*+K° "T °D+->D'OA-+ 



B D+ -£>•»*+ Jr D+-*D*o^+/DS ntag ' ^ C 




untae 



™D+ 1 ->D*+K° ^Df 1 ^D*+K°/D*+ ' J\ C ~ * D 



-cxtr 

Using f(c—>D* + ) and f(c—>D°) [30], recalculated with the updated values of the 

o 



branching fractions [1], and calculating the fragmentation fraction into untagged D 

= 39.8 ± 1.9(stat.) ± 1.5(syst.)t^(br.) %, 

where £>£ ) *+^£ ) o 7r + is the branching fraction of the decay D* + — > D°ix + (67.7% [1]) and 
the third uncertainty is due to the branching-fraction uncertainties, yields 

—^ = 2.3 ± 0.6(stat.) ± 0.3(syst.) 

in comparison with the world average value of 1.27 ± 0.21 [1]. Isospin invariance requires 
the matrix elements of the two measured Df ± decay modes to be the same, while an 
enhancement of the D*°K + final state is expected due to the larger phase space [41]. 

Assuming that the decay width of the D^ is saturated by the D*K final states, i.e. 



B n + ^n*+rso + o 



D+-*D*°K+ 



'si 



yields f(c — ► D^) (Table 3). The measured fragmentation fraction value agrees with 
those obtained in e + e~ annihilations and is above the prediction of the thermodynamical 
model [42]. 

The ratio for the two 1 + states 

f(c - D+)/f(c - D?) = 0.31 ± 0.06(stat.)ig:8t(syBt.) 

represents the strangeness-suppression factor for P-wave charm mesons. The measured 
value agrees with measurements of the strangeness-suppression factor for the lowest-mass 
charm mesons [13,30,47] and with the value of 0.3, used by default in simulations based 
on the Lund string fragmentation scheme [48]. 
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7 Search for the radially excited charm meson _D* /+ 

To search for the D*' + — > D* + 7t + tt~ decays, a D*' + candidate was formed by combin- 
ing each selected D* + candidate (Section 4.1) with two additional tracks with opposite 
charges. The additional tracks were assumed to be pions (ir^r), and were required to sat- 
isfy the pion dE/dx hypothesis with l n > 0.01 (Section 4). To reduce the combinatorial 
background, the cuts ^(tt^) < 1.1 and cos#*(.D* + ) < 0.8 were imposed, where 6*(D* + ) 
is the angle between the D* + in the D* + n^ir~ rest frame and the D* + 7r^n~ line of flight 
in the laboratory frame. To further reduce the combinatorial background, the following 
requirements were applied: 

Pt(^) > 0.15 GeV, p T (D* + 7r+7v;)/E^ >w > 0.25 

for the D* + decay channel (1) and 

p T (vr a ± ) > 0.25 GeV, p T (D*+7i+n~)/E e T >w > 0.30 

for the D* + decay channel channel (2). 

For each D*' + candidate, the extended mass difference, AM ext = M(Kttt: s tt"^7:~) — 
M(KiT7r s ) or AM ext = M{K'K'K'K'k s tx'^'k~) — M(Ktttttttt s ), was calculated. The invariant 
mass of the £>*+7r+7T~ system was calculated as M(£>* + 7T+7r~) = AM cxt + M(D* + ) PBG . 
The resolution in M(D* + 7i^7r~) around 2.64 GeV, where a narrow signal was reported by 
the DELPHI Collaboration [10], was estimated from MC simulations to be 5.6 MeV. 

Figure 10 shows the M(D* + 7r^n~) distribution below 2.9 GeV. The distribution was in- 
vestigated in the full accessible range; no narrow resonance was observed. 

An estimate of the fraction of D* + mesons originating from the D*' + — ■» D* + -k + tt~ de- 
cays was performed in the signal window of 2.59 < M(D* + 7i^n~) < 2.69 GeV. This 
window covers both theoretical predictions [9] and the DELPHI measurement [10]. The 
M(D* + 7t^ti~) distribution was fitted outside the signal window to the background func- 
tional form with two shape parameters, x A exp(—Bx), where x = AM ext — 2m n +. The 
number of reconstructed D*' + mesons was estimated to be 104 ± 83 by subtracting the 
background function, integrated over the signal window, from the observed number of 
candidates in the window. 

The number of reconstructed D*' + — > D* + 7t + tt~ decays was divided by the number of 
reconstructed D* + mesons, yielding a fraction of D* + mesons originating from the D*' + 
decays. To correct the measured fraction for detector effects, the relative acceptance was 
calculated using the MC simulation (Section 3) as a ratio of an acceptance for the D*' + — > 
D* + n + 7T~ state to the inclusive D* + acceptance. The acceptance of the requirement 
l„ > 0.01 for the additional tracks was calculated with data (Section 5.4). Subtraction 
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of the small 6-quark contribution, performed under a conservative assumption that all 
D*' + mesons are produced in charm fragmentation, changed the relative acceptance by 
~ 1.7% of its value. The relative acceptance was found to be 34% in the kinematic range 
described in Section 4.1. 

The fraction, JF, of D* + mesons originating from D*' + decays was calculated in the kine- 
matic range \rj(D* + )\ < 1.6 and pt{D* + ) > 1.35 GeV for the D* + decay channel (1), 
combined with channel (2) for pt{D* + ) > 2.8 GeV: 

Fd*>+^d*+^-/d*+ = 0.54 ± 0.43(stat.)i^(syst.) %. 

The fraction measured in the restricted pt{D* + ) and 7]{D* + ) kinematic range was ex- 
trapolated to the fraction in the full kinematic phase space (Section 5.4). Applying the 
estimated extrapolation factor, ~ 1.2, gives 

n x - + -,D* + ^-/D. + = °' 67 ± 0.53(stat.)+Ho(syst.) %. 
In the full kinematic phase space, the extrapolated ratio can be expressed as 

-r-cxtr _ J \ C ~^ ) U 

J ~D*'+^D*+tt+tt-/D*+ — J7~ — n*+N ' °D*'+^D*+tt+-k~-, 

where the fragmentation fraction /(c — > D*' + ) is the rate of c quarks hadronising as D*' + , 
and B£)*'+-^D*+n+n- is the branching fraction of the decay D*' + — > D* + tt + tt~ . 

Using /(c — > D* + ) [30], recalculated with the updated branching fractions [1], an upper 
limit was set on the product of the fraction of c quarks hadronising as a D*' + meson 
and the branching fraction of the D*' + — ► D* + 7i + tt~ decay in the mass range 2.59 < 
M(D*+7T+-K-) < 2.69 GeV: 

f{c -> D*'+) ■ B D , /+ ^ D , +w+n - < 0.4% (95% C.L.). 

The upper limit is the frequentist confidence bound calculated assuming a Gaussian prob- 
ability function in the unified approach [49]. It is stronger than the 0.9% limit on _D* /=t 
production in charm fragmentation obtained by OPAL [11]. 

The ratio of the D*' + -* D* + ix + ti- to £>?, Df -> L>* + 7r~ decay yields, calculated as 

-rextr 

„ _ J~ D*'+^D* + TT+TT-/D* + 

l<-D'"+->D*+7r+TT- /D 1 ,Df^D*+w- ~ :^xta , -r-cxtr ' 

-> D°->D*+n-/D*+ +^D*°^D*+7r-/D*+ 

is compared with those obtained by DELPHI [10] and OPAL [11] in Table 5. The 
ZEUS measurement is more sensitive to the existence of a narrow resonance decaying 
to D* + 7T + 7r~ . However, it is sensitive only to the resonance production in charm fragmen- 
tation while the LEP measurements are also sensitive to beauty fragmentation. 
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8 Systematic uncertainties 

The systematic uncertainties of the measured values were determined by varying the 
analysis procedure and repeating all calculations. The sizes of the variations were chosen 
commensurate with the estimated uncertainties of the relevant parameters and variables. 
The following groups of systematic uncertainties were considered. 

• {Si} The uncertainties related to the signal and helicity extraction procedures were 
obtained as follows: 

— for the D* + signals: the ranges for the background normalisation were reduced by 
2 MeV on either side; the fit was used instead of the subtraction procedure; 

— for the D + signal: the range for the signal fit was reduced by 20 MeV on either 
side; the amounts of the subtracted Df and A+ reflections were varied in the range 
of their uncertainties; a higher-order polynomial was included in the background 
parametrisation; 

— for the untagged D° signal: the range for the signal fit was reduced by 20 MeV 
on either side; the value of M(Ktt), where the background form with the expo- 
nential enhancement turns into the linear form, was varied between 1.84 GeV and 
1.88 GeV; a higher-order polynomial was included in the background parametrisa- 
tion; 

— for the D® and D^ signals: the ranges for the signal fit were reduced by 20 MeV 
on either side; higher-order polynomials were included in the exponential of the 
background parametrisations; the masses and widths of the wide excited charm 
mesons were varied in the range of their uncertainties [1] and their yields were 
varied by ±50%; 

— for the D\ helicity distribution: the acceptance dependence on the helicity angle 
was varied in the range of its uncertainty; the background functions in the four 
helicity intervals were allowed to have separate normalisations; 

— for the D^ signals: the ranges for the signal fit were reduced by 12 MeV on 
the upper side; higher-order polynomials were included in the exponential of the 
background parametrisations; the average shift of the signal in the M(D°K a ) 
distribution with respect to the mass of D^ x meson was varied in the range of its 
uncertainty (Section 6.2); 

— for the D^ helicity distribution: the acceptance dependence on the helicity angle 
was varied in the range of its uncertainty; the background function was allowed to 
have a free helicity parameter; 
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— for the D*' + signal search: the range for the background fit was reduced by 12 MeV 
on the upper side; a higher-order polynomial was included in the exponential of 
the background parametrisation; 

{#2} The uncertainty of the tracking reconstruction and simulation was taken into ac- 
count by varying all momenta by ±0.1% (magnetic field uncertainty) and by changing 
the track momentum and angular resolutions by ±5% of their values. 

{5 3 } The uncertainties of M(D*+) PDG - M(D+) PDG , M(D* + ) PDG - M(D°) PDG and 
M(K S ) PP)G were included. 

{$4} The uncertainties of the dE/dx requirements applied to the additional tracks 
(Sections 5.4, 6.4 and 7) were taken into account. 

{#5} The uncertainty of the CAL simulation was determined by varying the CAL 
energy scale by ±2%. 

{#6} The uncertainties of the fragmentation fractions f(c — > D* + ), f(c — > D + ) and 
/(c — > -Du n tag) were determined by adding in quadrature their statistical and system- 
atic uncertainties and the errors originating from the branching-fraction uncertainties. 
The uncertainty of the branching fraction of the K s decay into 7r + 7r~ [1] was also 
taken into account. 

{£7} The model dependence of the acceptance corrections was estimated by varying 
the p T (D* + , D + , D°) and r}{D* + , D + , D°) distributions of the MC sample by their 
uncertainties; the MC fraction of the lowest-mass charm mesons produced in a vector 
state was taken to be 0.6 ±0.1. 

{5$} The uncertainty of the beauty subtraction was determined by varying the b- 
quark cross section by a factor of two in the MC sample and by varying the branching 
fractions of 6-quarks to charm hadrons by their uncertainties [36-39] . 

{Sg} The extrapolation uncertainties were determined by varying relevant parameters 
of the Pythia simulation using the Bowler modification [27] of the Lund symmetric 
fragmentation function [28] 6 . The following variations were performed: 

— the mass of the c quark was taken to be 1.5 ± 0.2 GeV; 

— the strangeness suppression factor was taken to be 0.3 ± 0.1; 

— the fraction of the lowest-mass charm mesons produced in a vector state was taken 
to be 0.6 ±0.1; 
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An adequate use of the Peterson fragmentation function [50] for the extrapolation was not possible 
due to the absence of predictions or measurements of the Peterson parameter values for all involved 
charm mesons. Using the Peterson fragmentation function with the same parameter value (0.05) for 
all charm mesons increases the extrapolation factors by 10 — 25%. 
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— production rates of the excited charm and charm-strange mesons were varied by 
±50% around the central values tuned to reproduce the measured fractions of c 
quarks hadronising into D®, D^ or D^; 

— the Bowler fragmentation function parameter r c was varied from the predicted 
value 1 to 0.5; the a and b parameters of the Lund symmetric function were varied 
by ±20% around their default values [21]. 

Contributions from the different systematic uncertainties were calculated and added in 
quadrature separately for positive and negative variations. The results are given in Ta- 
bles 6-7. 

The relatively narrow AM, M(Kim) and M(Kn) ranges, used for the excited charm 
and charm-strange meson studies, selected only the central parts of the D* + , D + and 
D° signals, respectively (Section 4). It was checked that increasing the narrow ranges by 
25 — 50% produced no effect on the results beyond the expected statistical fluctuations. 
Similarly, no systematic shifts were found when removing the rj(7i a , K a ) < 1.1 requirement 
from the excited state selections (Sections 5.1, 5.2, 6.2 and 7). It was also checked that 
the D® width value cannot be significantly reduced by including an interference between 
the signal and background. 

9 Summary 

Sizeable production of the excited charm and charm-strange mesons was observed in ep 
interactions. The measured masses of the D®, D^ and D^ are in reasonable agreement 
with the world average values [1]. The measured D® width is 

r(D?) = 53.2 ± 7.2(stat.)+^(syst.) MeV 

which is above the world average value 20.4 ± 1.7 MeV [1]. 
The measured D® helicity parameter is 

h{Dl) = 5.9+f 7 (stat.)+^(syst.), 

which is inconsistent with the prediction of h = for a pure S'-wave decay of the 1 + state, 
and is consistent with the prediction of h = 3 for a pure D-wave decay. In the general 
case of D- and S'-wave mixing, the allowed region of the mixing parameters is consistent 
with the CLEO measurement [35] and marginally consistent with the BELLE result [4]. 

The measured Df x helicity parameter is 

h{D$) = -0.74±°f 7 (stat.)^(syst.). 
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This value is inconsistent with the prediction of h = 3 for a pure D-wave decay of the 1 + 
state, and is barely consistent with the prediction of h = for a pure S'-wave decay. The 
measurement suggests a significant contribution of both D- and S'-wave amplitudes to the 
.D 5 i(2536) + — > D* + Kg decay. The allowed region of the mixing parameters is consistent 
with the CLEO measurement [45] and with the BELLE result [46]. 

The ratios of the dominant D^ and D^ branching fractions are 



B 



D-f-^D+n_ = 2 g ± g( stat )+0.5( syst ^ 



^Df^D*+n 



B n+ — >n*0jc+ 

- 2.3 ± 0.6(stat.) ± 0.3(syst.) 



>D+ 1 ^D*+K° 



in agreement with the world average values [1]. 

The fractions of c quarks hadronising into D®, D^ or D^ mesons are consistent with 
those obtained in e + e~ annihilations (Table 3), in agreement with charm fragmentation 
universality. Sizeable fractions of the D* + , D + and D° mesons emanate from these excited 
states. 

No radially excited D*' + meson was observed. An upper limit, stronger than that obtained 
by OPAL [11], was set on the product of the fraction of c quarks hadronising as a D*' + 
meson and the branching fraction of the D* /+ — > D* + tt + ti~ decay in the range of the D*' + 
mass from 2.59 to 2.69 GeV: 

f(c -> D*'+) ■ B D * I+ ^ D *+^- < 0.4% (95% C.L.). 



Appendix 

Relativistic Breit-Wigner function 

The mass distribution, M, of a resonance with a non-negligible natural width decaying into 

two particles is described by a relativistic Breit-Wigner function with a mass-dependent 

width [51]: 

dN MM T(M) 

DC 

dM (M 2 - M 2 ) 2 + M 2 r 2 (M) ' 
M fp*\ 2l+1 .,, .. 



where Tq is the nominal resonance width, p* is the momentum of the decay products in 
the resonance rest frame and Pq is the value of p* at the resonance nominal mass M . 
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The hadron transition form- factor, F l (p*,p^), in the Blatt-Weisskopf parametrisation [52] 
equals 1 for S-wave (/ = 0) decays and 

KP ,Po) 9 + 3(p* r )2 + (p* r )4 

for D-w&ve (I = 2) decays, where r = 1.6 GeV -1 is a hadron scale. 
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decay 


£>* + channel (1) 


D* + channel (2) 


p T (K) (GeV) 


>0.45 


>0.5 


p T (n) (GeV) 


>0.45 


>0.2 


p T (ir s ) (GeV) 


>0.1 


>0.15 


p T (D*+)/E° T > 10 ° 


>0.12 


>0.2 


p T {D* + ) (GeV) 


> 1.35 


>2.8 


\ V (D*+)\ 


: 1.6 


: 1.6 


M(D°) (GeV) for 
p T (D* + ) < 3.25 GeV 


1.83-1.90 


1.845 - 1.885 


M(D°) (GeV) for 
3.25 < p T (D*+) <5GeV 


1.82-1.91 


1.845 - 1.885 


M(D°) (GeV) for 
5 < p T (D* + ) < 8 GeV 


1.81-1.92 


1.835 - 1.895 


M(D°) (GeV) for 
p T (D* + ) > 8 GeV 


1.80-1.93 


1.825-1.905 



Table 1: Requirements applied for selections of D* + candidates in the decay chan- 
nels (1) and (2) (see text). The mass resolution dependence on pt{D* + ) is taken 
into account in the requirement on consistency of the reconstructed and nominal 
D° masses. 
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final state 


D* + n a 


D + n a 


Signal yields 


NW) 


3110 ±340 




N(D*°) 


870 ± 170 


690 ± 160 


Background parameters 


Yield 


169 ±18 


1540 ± 300 


A 


0.37 ±0.3 


1.27 ±0.7 


B 


1.3 ±0.3 


7.7 ±0.4 


C 


-1.4 ±0.3 


2.3 ±0.3 



Table 2: The numbers of reconstructed D® and D^ mesons and values of all free 
background parameters yielded by the simultaneous fit of the M(D + ir a ) distribution 
and the M(D* + ir a ) distributions in four helicity intervals (see text). The mass, 
width and helicity parameters are given in the text. 





f(c - - £>?) [%} 


f(c - ♦ Df) [%] 


/(c - * -Di) [%] 


ZEUS 


3.5 ±0.4!*$ 


3.8 ±0.7+^6 


1.11 ±0.16lg;?g 


OPAL [38] 


2.1±0.7±0.3 


5.2 ±2.2 ±1.3 


1.6±0.4±0.3 


ALEPH [39] 






0.94 ±0.22 ±0.07 


Model [42] 


1.7 


2.4 


0.54 



Table 3: The fractions of c quarks hadronising into the D®, D^ and D~£± mesons 
(Sections 5.4 and 6.4)- The first uncertainty is statistical and the second is sys- 
tematic (Section 8). 



final state 


D* + K° s 


D°K a 


Signal yields 


N(D+) 


100 ±13 


136 ±27 


Background parameters 


A 


0.43 ±0.06 


0.43 ±0.05 


B 




4.3 ±1.0 



Table 4: The numbers of reconstructed Df r mesons and values of all free back- 
ground parameters yielded by the unbinned likelihood fit performed simultaneously 
using values of M(D°K a ), M(D* + K° S ) and helicity angle for D* + Kg combinations 
(see text). The mass, width and helicity parameters are given in the text. 
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T^D*'+^D*+tt+tt- /Dl,Df^D*+TT- 


DELPHI [10], Z° -^bb,cc 


49 ± 18 ± 10 % 


OPAL [11], Z° ^bb,cc 


5 ±10 ±0.2% 
<22% (95% C.L.) 


ZEUS, ep -> ccX 


4.5 ±3.6+°;?% 
< 12% (95% C.L.) 



Table 5: The ratio of the D*'+ -»■ £>* + 7T + 7r and D?, D|° -► £>* + 7T decay yzetas, 
^•D*'+-+D*+7r+7r-/D ,D* -»o*+7r-- 27* e /^ rs ^ uncertainty is statistical and the second is 
systematic (Section 8). 





total 


Si 


S 2 


S3 


M(D°) [MeV] 


±0.9 


+0.4 
-0.5 


±0.8 


±0.0 


M(Df) [MeV] 


+1.2 
-1.3 


+0.6 
-0.8 


±1.0 


+0.1 
-0.0 


r( J D°) [MeV] 


+3.3 
-4.9 


+3.3 
-4.9 


±0.2 


±0.0 


h(D<{) 


+2.4 
-1.0 


+2.4 
-1.0 


±0.0 


±0.0 


M(D+) [MeV] 


±0.10 


+0.06 
-0.05 


±0.08 


±0.02 


M^ + i) 


+0.06 
-0.05 


+0.06 
-0.05 


— 


±0.00 



Table 6: The total and S1-S3 (see text) systematic uncertainties for the mass, 
width and helicity parameters of the excited charm and charm-strange mesons. 
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total 

(%) 


6i 

(%) 


52 


(%) 


6, 
(%) 


s 5 

(%) 


^6 

(%) 


5 7 


s 8 

(%) 


69 
(%) 


•rcxtr 

X D°->D*+n-/D*+ 


+9.3 
-14.4 


+8.5 
-13.9 


+0.6 
-0.3 


±0.0 


±0.1 


+2.2 
-2.3 


— 


+1.1 

-0.6 


±0.7 


±2.6 


J~ Df->D*+w-/D*+ 


+6.5 
-7.1 


+5.1 
-5.9 


+0.3 
-0.5 


±0.0 


±0.1 


+2.4 
-2.1 


— 


+1.0 
-0.6 


±1.2 


+2.8 
-2.9 


-rextr 

J ~D*°->D+tt-/D+ 


+12.3 
-16.7 


+10.8 
-15.8 


+3.0 
-0.7 


+0.2 
-1.0 


±0.1 


+2.8 
-3.1 


— 


+1.0 
-0.4 


+1.4 
-1.0 


+4.6 
-4.2 


^Df^D+TT- 


+18.3 
-20.0 


+12.0 
-16.1 


+1.7 
-0.4 


±0.2 


±0.0 


+0.4 
-1.0 


+13.2 
-11.2 


+0.5 
-0.8 


+1.3 
-0.9 


+3.2 
-4.8 


"z>*0 _►£>*+*■- 


/(c-I>?) 


+11.5 
-16.4 


+8.5 
-13.9 


+0.6 
-0.3 


±0.0 


±0.1 


+2.2 
-2.3 


+6.9 

-7.8 


+1.1 
-0.6 


±0.7 


±2.6 


/(c - ^°) 


+12.3 
-14.6 


+ 8.2 
-11.8 


+0.9 
-0.0 


+0.1 
-0.7 


±0.1 


+2.7 
-2.8 


+7.7 
-7.1 


+0.3 
-0.0 


+1.2 
-1.0 


+4.0 
-3.6 


J ~Df 1 ^D*+K°/D*+ 


+4.5 
-4.1 


+1.6 
-2.0 


+0.7 
-0.3 


±0.0 


±0.0 


+0.1 
-0.0 


±0.1 


+1.7 
-1.0 


±0.6 


+3.7 
-3.3 


•r-cxtr 

D + — ►n*OR _ + /D° 


+6.3 
-8.3 


+1.9 
-4.0 


+3.0 
-0.7 


±0.2 


±0.3 


+3.8 
-3.5 


— 


+0.5 
-0.4 


+1.4 
-0.7 


+3.3 
-6.2 


ft 

si 


+12.5 
-13.3 


+2.7 
-4.3 


+2.6 
-0.4 


+0.2 
-1.0 


±0.3 


+3.6 
-3.5 


+ 11.0 
-10.3 


+1.2 
-1.8 


+0.9 
-0.4 


+2.0 
-6.0 


U D + -> D *+ K 

si 


f(c - D+) 


+7.4 
-8.6 


+1.4 
-2.8 


+2.3 
-0.6 


+0.1 
-0.2 


±0.2 


+2.6 
-2.8 


+5.3 

-6.1 


+0.6 
-0.4 


+1.1 
-0.7 


+3.2 
-4.8 


/r-extr 

J~ D*'+^D* + TT+TT-/D* + 


+4.7 
-15.0 


+3.0 
-13.9 


+1.6 
-2.0 


— 


±0.2 


+2.0 
-2.4 


— 


+1.3 

-0.7 


+2.0 
-1.0 


+1.2 
-4.5 



Table 7: The total and Si-Sg (see text) systematic uncertainties for extrapolated 
fractions, for ratios of the dominant branching fractions and for fragmentation 
fractions of the excited charm and charm-strange mesons. 
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Figure 1: The distributions of the mass differences (dots), (a) AM = M(K7rir s ) — 
M(Ktt) for D*± -> (K-n)-K s candidates and (b) AM = M{Kmx-mi s ) - M(Kixixix) 
for D*^ — > (A'7T7T7r)7r s candidates. The solid curves represent fits to the sum of a 
modified Gaussian function and a background function. The histograms show the 
AM distributions for wrong-charge combinations. Only D* ± candidates from the 
shaded ranges were used for the analysis of excited states. 
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Figure 2: The distributions of the invariant masses (dots) for (a) the D^ — > K7nr 
candidates and (b) the D°/D° — > Kir candidates after the reflection subtractions. 
The solid curves represent fits to the sum of a modified Gaussian function and a 
background function (dashed curves). Only candidates from the shaded ranges were 
used for the analysis of excited states. 
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Figure 3: The distribution of (a) M(L>* ± vr a ) = AM cxt + M(D* + ) PDG , where 
AM ext = M{Kinr s 7T a ) - M{Ktttt s ) or AM cxt = M(Kimivn s ii a ) - M{K7nnnr s ), 
for £>?, D*° -► £>*% candidates and (b) M{D ± Ti a ) = AM cxt + M( J D+) PDG , where 
AM ext = M(K7T7T7r) - M(K-K7r), for Df -> £>% candidates (dots). The solid 
curves represent the result of the simultaneous fit with the background contribution 
given by the dashed curves (Section 5.3). Contributions from the wide -Di(2430)° 
and -Dq(2400)° states are shown in (a) and (b), respectively, as shaded bands. The 
histograms show the distributions for wrong-charge combinations. 
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Figure 4: The distribution of ' M{D* ± n ) = AM cxt + M(D* + ) PDG for D\, D 
D^it candidates in four helicity intervals: (a) \ cosa| < 0.25, (b) 0.25 < | cosa| < 
0.5, (c) 0.5 < |cosa| < 0.75 and (d) |cosa| > 0.75 (dots). The solid curves 
represent the result of the simultaneous fit with the background contribution given 
by the dashed curves (see text). 
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Figure 5: Cosine of the relative phase of S- and D-wave amplitudes versus 
r = T S /{T S + V D ) in the Di(2420)° -> D* + n~ decay from the ZEUS, CLEO and 
BELLE measurements. There is a marginal overlap between the ranges defined by 
the ZEUS and CLEO measurements. The difference between the ZEUS and BELLE 
measurements, evaluated with Eq. (6), is ~ 2 standard deviations. 
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Figure 6: The distribution of the invariant mass, M(it + ii~), in events with a 
D*^ candidate. The solid curve represents a fit to the sum of a modified Gaussian 
function and a linear background function (dashed curve). Only Kg candidates from 
the shaded range were used for the analysis of the excited charm-strange mesons. 



39 



ZEUS 



O 80 
^ 70 



10 




a 



• ZEUS (126 pb" 1 ) 

— simultaneous fit 
--- background fit 

N(D^) = 100±13 




t 



T 



T 



1 

1 

(1 

1 


1 
1 

1 


' 1 ' 


n 


i ii 


ii 

i 



2.52 2.54 2.56 2.58 2.6 2.62 2.64 

M(D*K9 = AM ext + M(D* + ) PDfi + M(K°) PDR (GeV) 



N(D* 1 ) = 136 + 27 




M(D°K a ) = AM ext + M(D°) PDG (GeV) 

Figure 7: The distribution of (a) M(D* ± iT§) = AM ext + M(D*+) PDG + 
M{K S ) PDG , where AM cxt = M{K7nr s 7r + 7T-) - M{Ktt7t s ) - M(tt + tt~) or AM cxt = 
M(K nmm s n + 7i~) — M(K7nr7nr s ) — M(ti + 7i~), for Df x — > D* ± K° S candidates and 
(b) M{D°K a ) =_AM cxt + M(D°) PDG , where AM cxt = M(KnK a ) - M(Ktt), for 
D sl — >• D*°K + /D*°K~ candidates (dots). The solid curves represent the result of 
the simultaneous fit with the background contribution given by the dashed curves 
(Section 6.3). 
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Figure 8: The distribution of M(D m± K^) = AM ext + M(D* + ) PDG + M(K° S ) PDG 
for Df x — > D* ± K < g candidates in four helicity intervals: (a) |cosa| < 0.25, (b) 
0.25 < |cosa;| < 0.5, (c) 0.5 < |cosa| < 0.75 and (d) |cosa| > 0.75 (dots). 
The solid curves represent the result of the simultaneous fit with the background 
contribution given by the dashed curves (see text). 
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Figure 9: Cosine of the relative phase of S- and D-wave amplitudes versus 
r = T S /(T S + T D ) in the D sl (2536)+ -> D*+K° s decay from the ZEUS, CLEO and 
BELLE measurements. 
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Figure 10: The distribution of M(D* ± 7r+7r~) = AM cxt + M(D*+) PDG , 

where AM ext = M(Kir7r s ir+7T~) - M(Kinr s ) or AM ext = M(K'innnr s ir+'K-) - 
M(K7nnnr s ), for D* r± —> D* ±r K +r K~ candidates (dots). The inset shows the D* f± 
signal window covering both theoretical predictions and the DELPHI measurement. 
The solid curve is a fit to the background function outside the signal window. The 
shaded histogram shows the Monte Carlo D* l:k signal, normalised to the obtained 
upper limit (95% C.L.) and shown on top of the fit interpolation (dashed curve). 
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